Surface Contact Potential Patches and Casimir Force Measurements 
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We present calculations of contact potential surface patch effects that simplify previous treat- 
ments. It is shown that, because of the linearity of Laplace's equation, the presence of patch po- 
tentials does not affect an electrostatic calibration of a two-plate Casimir measurement apparatus. 
Using models that include long-range variations in the contact potential across the plate surfaces, a 
number of experimental observations can be reproduced and explained. For these models, numerical 
calculations show that if a voltage is applied between the plates which minimizes the force, a resid- 
ual electrostatic force persists, and that the minimizing potential varies with distance. The residual 
force can be described by a fit to a simple two-parameter function involving the minimizing potential 
and its variation with distance. We show the origin of this residual force by use of a simple parallel 
capacitor model. Finally, the implications of a residual force that varies in a manner different from 
1/d on the accuracy of previous Casimir measurements is discussed. 

PACS numbers: 31.30.jh 12.20.-m 42.50.Ct, 12.20.-m, 78.20.Ci 



I. INTRODUCTION 

It is often assumed that the surface of a conductor is 
an equipotential. While this would be true for a perfectly 
clean surface of a homogeneous conductor cut along one 
of its crystalline planes, it is not the case for any real 
surface. Potential patches can be caused by, for example, 
oxide films or some other films adsorbed on the surface 
(less than a monolayer is required), strains in the sur- 
face, and chemical impurities within the surface. Such 
patches effectively create a surface dipole layer, which 
alters the potential above the surface. Various types of 
monopolar charge disorder could be also present, signifi- 
cantly altering the properties of a surface from that of an 
ideal conductor [l[. Even for chemically unreactive no- 
ble metals, such as gold and copper, carefully prepared 
in an ultra-clean environment to minimize such "dirt" 
films, experiments show that typical surface potential 
variations are on the order of at least a few millivolts 0- 
[ij]. Similar effects were found in a more recent measure- 
ment on a pair of metallic plates employed in the Laser 
Interferometer Gravitational Wave Observatory (LIGO) 
project The cause of surface potential variations is 
most likely local changes in surface crystalline structure, 
giving rise to varying work functions and hence varying- 
potcntial patches. It is well known that the work function 
of a metal surface depends on the crystallographic plane 
along which it lies; as an example, for gold the work func- 
tions are 5.47 eV, 5.37 eV, and 5.31 eV for surfaces in the 
(100), (110), and (111) directions, respectively @. This 
variation is most likely due to the different effective elec- 
tron masses, and the resulting different Fermi energies 
along the corresponding directions. 

When two conductors of different work functions are 
brought into contact, electrons flow until the chemical po- 
tential (i.e., the Fermi energy) in both conductors equal- 
izes. As a result, there is a net dipole distribution created 
at the interface, which gives rise to a "contact" or "Volta" 



potential established between the two conductors, equal 
to the difference between their work functions. We there- 
fore expect that the electrostatic potential along a chem- 
ically clean metal surface varies on the length scale of the 
typical size of surface crystallites, which can vary from 
sub-micron to millimeter or larger scales. 



The generation of an attractive force between conduct- 
ing surfaces due to metallic contacts has long been rec- 
ognized as a possible systematic limitation to a Casimir 
force measurement 0, Q ■ Recently, Spcake and Trenkel 
have performed a formal treatment of the effect of 
random, zero-average, patch potentials. We reproduce 
here their result using a slightly different formalism, and 
show that the result can be simply obtained, expressed 
in terms of the surface potential autocorrelation func- 
tion. Furthermore, we provide a simple model that shows 
the origin of a distance-dependent minimizing potential, 
and that this dependence leads to a residual electrostatic 
force that can have a complicated character. The goal of 
this paper is not to present a full rigorous mathematical 
derivation of these effects, but to illustrate their funda- 
mental nature and motivate functional forms that have 
been observed in experiments in a straightforward man- 
ner. More importantly, if the patch potential is actually 
measured, it would be possible to perform an exact nu- 
merical calculation of the excess force. Our recent work 
[Tol | indicates variations of contact potentials on the level 
of 10 mV, which is a challenging level for Kelvin probe 
techniques with the appropriate spatial resolution fill ]. 
Thus, the full experimental description of these effects 
remains an open challenge. 



In this paper we define a positive force as an attrac- 
tion between the surfaces under consideration. There- 
fore, attractive electrostatic and Casimir forces are taken 
as positive. 
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II. ATTRACTIVE ELECTROSTATIC FORCE 
DUE TO SMALL SURFACE PATCHES 

In this section we will consider the effect of random 
surface patches on the electrostatic interaction between 
metallic plates, for the plane-plane and plane-sphere ge- 
ometries. We will assume that the typical patch area is 
much smaller than the effective area of the surface inter- 
actions, defined for each of the geometries below. 



A. Parallel planes geometry 

Consider two plane parallel metallic surfaces at z = 
and z = d. Let the electrostatic potential at z = be 
V(x,y,z = 0) = V a (x, y) and at z = d be V(x,y, z = 
d) = Vb{x 1 y). The solution of Laplace's equation in the 
region < z < d can be easily found separating variables 
in cartesian coordinates, V(x, y,z) — X (x)Y (y) Z (z) , 
where 



1 d 2 X 
X dx 2 



Y~dtf 



1 d 2 Z 
'Z~dz l: 



= 7 2 , 



(1) 



and 7 2 = a 2 + fP . The general solution of Laplace's 
equation in this geometry can be written as 



In order to proceed, we will assume that the electrostatic 
patches are stochastic, uncorrelated between the differ- 
ent plates, and for a given plane the 2-point correlation 
function is diagonal in the k-basis. That is 



(K.k) - (H,k> = (%,kK,k-) = 0; 
(K,kK,k') = C a , k <5 2 (k-k'); 
(H.kH.k/) = C feik <5 2 (k-k'), 



(6) 



where (...) means stochastic average. Note that, as in [§], 
we assume zero cross-correlation between the plates, and 
that the correlation function of each plate is independent 
of the position of the other plate (i.e., independent of the 
distance d). 

It is then easy to calculate each of the terms in the 
expression for U pp : 



o > 



d 2 k k 2 sin 2 (k x x) cos 2 (k y y) 



(2tt) 2 4sinh 2 ( 7 d) 
x {2CVk[cosh(2 7 z) - 1] + 2C a , k [cosh(2 7 (z - d)) - 1]} 



dV 
dy 



) 



d 2 k k 2 cos 2 (k x x) sin 2 (k y y) 



(2tt) 2 4sinh 2 ( 7 d) 
x {2C 6jk [cosh(2 7 z) - 1] + 2C a , k [cosh(2 7 (z - d)) 



1]} 



V (x, y,z) = / da dp [cos(ax) + A a sin(ax)] 

J —OQ 

x [cosC%) + Ap sinC%)] [B+e** + e^ z ] . (2) 

It is convenient to expand the boundary conditions in a 
cosine Fourier series: 



V a (x,y) 



d 2 k 



VaMCOs(k x x)cos(k y y), (3) 



and similarly for V b (x, y). Imposing the boundary con- 
ditions we get A a = Ap = 0, a = k x , f3 = k y (hence 

7 = y^Tfc 2 ), B+ = (V fc , 7 - K, 7 e-^ d )/2sinh( 7 d) and 

B~ = (V^. 7 e 7d — V&. 7 )/2 sinh( 7 <i). Therefore, the general 
solution of Laplace's equation with the given boundary 
conditions on each plate is 



V(x,y,z) = 
x \e" (V bJ 



d 2 k cos(k x x) cos(k y y) 
(2tt) 2 2 sinh( 7 d) 
K, k e"^) + e-T* (K, k e 7d - V blk )] 



(4) 



Now we calculate the electrostatic energy between the 
plates. The electrostatic energy density is u = y|E| 2 , so 
the total energy U pp per unit area A = L x L y is 



pp 
2 A 



dx 



dV 

dx 



L y rd 

dy I dz 
Jo 
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dV 
dy 



dV 

dz 



(5) 



> 



d 2 k 7 2 cos 2 (fc x x) cos 2 (k y y) 



((-, 

\dz J ' J (2tt) 2 4sinh"( 7 d) 
x {2C 6 , k [cosh(2 7 z) + 1] + 2C a , k [cosh(2 7 (z - d)) + 1]} 



The 



and 



y 



integrations are trivial, since 
L x 1 Jq L * dx sin 2 (fc x a;) = L^ 1 J^ x dx cos 2 (k x x) = 1/2 
and the same for y. For the z integration we 
use Jq Gk[cosh(2 7 z) ± 1] = sinh(2 7 d)/2 7 ± d and 

J d dz[cosh(2 7 (z - d)) ± 1] = sinh(2 7 d)/2 7 ± d. Finally 
the electrostatic parallel-plate energy is 



U 

pp ~ 16 



d 2 k 7 sinh(2 7 d) 
(2tt) 2 sinh 2 ( 7 d) 



[Ca, 



b.kj 



(7) 



In the special case of an isotropic patch distribution, the 
correlation functions depend only on k = |k|, that is 
C a .k = C a h and C b u = C b . k - Then 



U, 



16 27r7 sinh 2 (fcd) L ' 



Cb,k 



(8) 



For the plane-plane geometry, the effective area of in- 
teraction A e g is the same as the total area of the plate 
A c ff = A. In this case, the small surface patch limit 
corresponds to k 2 A 1, which basically means that we 
neglect finite-size effects in the computation of the elec- 
trostatic energy. 

Let us analyze different limiting cases of Eq.©. For 
surface potential patches small with respect to the effec- 
tive interaction area (k 2 A 3> 1) but large with respect 
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to the plates separation (kd <C 1 or d <c A, where A is 
a characteristic length of a potential patch) the energy 
scales as 1/d, which is the same as in the usual electro- 
static case for fixed (non-stochastic) potential difference. 
This limit kd — > is essentially the "proximity force ap- 
proximation" (PFA), to be discussed more fully in the 
next Section, applied to the electrostatic problem with 
in-plane potential variations. In this approximation the 
net energy and force is calculated by considering the at- 
traction between paired infinitesimal surface elements on 
each plate, and for large surface patches, the energy is 
then the additive sum of the usual 1/d plane-plane en- 
ergies. Indeed, if we rewrite Eq. (8) defining the root 
mean square (rms) potential fluctuations V 2 lls as 

1 poo poo 

KL = / dk k(C a . k + C b , k ) = / dk kS(k), (9) 
on Jo Jo 

one obtains in the limit kd — > 0, 

eoV 2 

u PP = -^p. (10) 

In the opposite limit (surface potential patches small 
with respect to the effective interaction area, k 2 A ^> 1, 
and small with respect to the plates separation, kd ^> 
1) Eq.(j8]) has an asymptotic behavior independent of 
the distance d. This is an artifact of the calculational 
method, that has included the self-energy of each plate. 
Following we remove from the above the potential 
energy at infinite separation, in order to have an ex- 
pression for the interaction energy only. Using that 
sinh(2fcd) = 2sxah(kd) cosh(fcd), we get 

fn f°° U 2 P~ kd 

= 17 / dk f- - S(k). (11) 
2 J sinh(fcd) w v ' 

Therefore, in the limit kd — >• oo the interaction energy 
vanishes exponentially. The reason is that in this case 
the patches are so small and change sign so rapidly that 
there is no net electric field at a large distance d from a 
given plate, hence there is no interaction with the other 
plate. 

Some remarks are in order. First, we have assumed 
that the boundary conditions on the two plates had only 
stochastic components fluctuating around 0. When there 
is an external fixed potential difference V between the 
two plates, the energy is the sum of the usual V 2 term 
plus the patch component calculated above (basically this 
is due to the linearity of Laplace's equation). The cross- 
terms VV^k and VV^k cancel upon taking stochastic av- 
erage. Second, when the cross-correlations (Vi^V^k') are 
not zero, there is a possibility that the interaction energy 
depends on the relative lateral position of the two plates, 
and hence it is possible to have a lateral force between 
the plates due to stochastic patches. 



B. Sphere-plane geometry 

In order to compute the patch effect on the force in the 
sphere-plane configuration we make use of the proximity 
force approximation. Just as in the case of roughness in 
Casimir physics, one must distinguish between two PFAs: 
one is for the treatment of the curvature of the sphere 
(valid when d R, where R is the radius of curvature), 
and the other one is the PFA applied to the surface patch 
distribution (valid when kd <C 1). We assume that we 
are in the conditions for PFA for the curvature, but we 
keep kd arbitrary. 

For the sphere-plane geometry we define an effective 
area of interaction by calculating along the plane surface 
the distance r from the point of closest approach (r = 0, 
corresponding to the minimal sphere-plane separation d) 
where the separation between the surfaces doubles. For 
a given r, this latter separation is given by d(r) = d + 
fl(l-cos#), with sin 6 = r/R. In the limit 9 < 1 (which 
corresponds to R ~3> d(r)) we obtain d(r) = d + r 2 /2R. 
The condition for the surface separation to be double 
that of the closest-approach distance is d(r e s) = 2d = 
d + r 2 s /2R, and the effective area is then 

A cS = nr 2 B = 2nRd. (12) 

For a sphere of radius R = 15 cm separated by a plane 
by a distance d = 1 /im, the effective distance above 
defined is r e g = 0.05 cm and the effective area is A e fi = 
0.009 cm 2 . The small patch limit corresponds to surface 
patches of area much smaller than this effective area A e g , 
that is k 2 A cS > 1. 

In the proximity force approximation, the electrostatic 
force in the sphere-plane case is F sp (d) = 2tt RU pp (d), 
namely 

poo >2 -kd 

F sp = neoR / dk——S(k). (13) 
Jo sinh(fea) 

There are a number of models that can be used to de- 
scribe the surface fluctuations. The simplest is to say 
that the potential autocorrelation function is, for a dis- 
tance r along a plate surface, 

K(r) = V^e.- r2 / x ' 2 . (14) 

Then, by the Wiener-Khinchin theorem, the power spec- 
tral density S(k) can be evaluated as the cosine two- 
dimensional Fourier transform of the autocorrelation 
function, which in our notation is [12J 

S{k) = 2-V 2 ^e~^ k \ (15) 

where the factor of two reflects the statistically- 
independent contributions from the two plate surfaces 
(Ci,fc + C2,fc)- The plane-sphere force is then given by, 
using k = u/X, 

y2 poo e _7r " 2 

F sp = 2ne R^- duu \ 2nd/x - ■ (16) 



4 



The limit of small potential patches kd — > (and also 
small with respect to the effective area of interaction, 
k 2 A e g ;§>!), the force is 



e a RV 2 



(17) 



suggesting that V 2 ms = V 2 /it. For the large kd limit, the 
force becomes exponentially small, just as in the parallel- 
plates geometry. 

Another possible model for the patch distribution is 
the one used in Assuming that C a ,k = CVfc = V 2 = 
const for k nlln < k < /c max and zero otherwise, we get an 
expression for V 2 ms (similar to Eq. (12) in |[) 



V 



1 



V 2 



— dkk(C a<k + C btk ) = ^-(k 



dkkS(k). 



We then obtain 



S(k) = 



V 2 
1 o 



2K 



2 

rms 



4-tt k' 2 — k 2 ' 
^" ""max "'min 



k 2 ) 



(18) 



(19) 



for fc m in < k < fc max and zero otherwise. The sphere- 
plane force from Eq. ([15]) therefore is 



F. 



4irenV 2 R 



S P 1,2 _ 1,2 

'"max %iin J k„ 



2 c -kd 



dk- 



k 2 e 



sinh(fcd) ' 



(20) 



which is the identical (apart from an overall, conventional 
sign) to Eq. (14) of jg. 



III. ATTRACTIVE FORCE DUE TO LARGE 
SURFACE PATCHES 

When the surface patches are larger than the effective 
area of interaction A c r defined in the previous section, 
the force between the plane and the sphere due to elec- 
trostatic patches can still be calculated using the method 
described above, but the average potential, as inferred by 
measuring the voltage at which a minimum in the attrac- 
tive electrostatic force occurs, will vary with distance. 
This can be thought of as a finite-size effect; if the patch 
size is roughly the diameter of the plates, then there can 
be a non-zero average over the surface. Alternatively, this 
problem can be addressed by assuming a slowly varying 
average potential across the plate surfaces, as developed 
below. 

In the PEA, the plate surfaces are divided into differ- 
ential areas, and the attractive force, given the potential 
difference between the plates and the derivative of the 
capacitance between them, is calculated. Specifically, for 
the sphere-plate geometry, 



r 2ir 



dip 



rdr 



(d + r 2 /2R) 21 



(21) 



where V(r, (f) is the net potential difference between the 
surface differentials, located at (r, ip) relative to the point 
of closest approach (r = 0). Writing V in this form allows 
the possibility that there can be a slow (coherent) varia- 
tion across the surface, as opposed to, and in addition to, 
innumerable small random patches. The PEA works here 
for two reasons. One is that the lines of electric force do 
not cross each other, the second is that the radius of cur- 
vature is large, so the angular deviations of the field lines 
are small when the plate diameter D satisfies D <C R- In 
this limit, the assumption that each surface differential 
element interacts only with a single element in the other 
plate is a good approximation. 

A slow variation in potential across the plate surfaces 
will manifest itself as a distance variation in the poten- 
tial that minimizes the electrostatic attractive force, i.e., 
V m = V m (d). Specifically, if we define the force with 
some externally applied voltage Vq to be 



p2ir 



Va 



I rar ( d + r 2 /2R) 2 



(22) 



the minimized force at a fixed distance determines the 
minimizing potential, 







dF(d,V ) 



dV 



V =V m 



2tt 



d<p 



rdr 



V{r, ip) + V m 
{d + r 2 /2Rf 



This equation implies a minimizing potential dependent 
on distance, V m = V m (d). Note that, in the idealized 
case of an equipotential surface, i.e., V(r, ip) = const, 
V m would be independent of d, and the minimized elec- 
trostatic force F(d, Vq = V m ) vanishes. Incidentally, the 
second derivative of F(d, Vq) with respect to Vb can be 
used to determine the distance at which the measurement 
is being made, 



d 2 F(d, Vq) 
dV 2 



= 2tt£ 



dr 



2ttR6, 



(d + r 2 /2R) 2 



where the finite size effects are neglected (upper limit of 
the r integration is set to infinity, which is a very good 
approximation when d <C R). The important implication 
is that the patch potentials do not interfere with the elec- 
trostatic calibration, that is the fundamental basis of our 
experiment (Toj . and of all Casimir force experiments. 

It is worth emphasizing a couple of points. First, that 
the origin of the distance dependence of the minimizing 
potential V m (d) is an interplay between the curvature of 
the surfaces and a variation of the electrostatic potentials 
Vi (i = a, b) along the surfaces (possibly due large surface 
patches). We have shown above how this effect arises in 
the context of the sphere-plane geometry, but of course 
it can be easily generalized to any geometry involving 
non-planar surfaces. It also follows from the above that 
for the parallel plates geometry one should expect that, 
even when the surface potentials Vi = Vi(x, y) vary along 
the plane surfaces, there should be no interplay with the 
(infinite) curvature of the planes, and thus the minimiz- 
ing potential V m should be distance-independent [l3j]. 
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G 



t t 



Cn 



implying a residual electrostatic force 



F T t(d) = F{d,V = V n , 

(-112 
C 4- — — 



C' + C' 



2 

V C 

2 ' 



(26) 



FIG. 1: A toy model illustrating the mechanism for the gen- 
eration of a distance-dependent minimizing electrostatic po- 
tential Vm(d) and electrostatic residual force F^ s (d). 



It is easy to take a case of parallel plate capacitors (C' a = 
—eoA/d 2 and C' b = —e A/ (d+ A) 2 , where A is the area of 
each of the upper plates in Fig. 1, assumed to be equal; 
hence, the lower continuous plate has 2A) and to 
show that there is a residual electrostatic force at the 
minimizing potential. Indeed, in such case, 



Second, distance dependence of the electrical potential 
minimizing the force between the plates has been ob- 
served in a number of experiments in the sphere-plane 
geometry [l4l - fl6j . as well as in our own work [loj . with 
further investigations under way. 

It has been suggested that the variation of the min- 
imizing potential with distance can cause an additional 
electrostatic force F(d, Vb = V m (d)), and an estimate 
was made for the possible size of the effect [l7|, where 
the varying contact potential is considered in a system of 
plates connected in series. The analysis presented in [ijj 
does not reproduce the effects seen in our experimen- 
tal work, and we were unable to develop a fundamen- 
tal theory of a plate-plate interaction that could cause 
a varying contact potential. Nonetheless, we have seen 
above that within the PEA a coherent variation of the 
surface potential along the non-planar surface does im- 
ply a distance-dependent minimizing potential. Now we 
present a simple model that produces not only varying 
contact potentials, but also the corresponding residual 
electrostatic force, consistent with our observations in 
[lOl] . The model is depicted in Fig. [1] In this figure, 
the two capacitors (short distance, C a (d), long distance, 
Cb{d + A)) create a net force on the lower continuous 
plate (setting Vi = initially), 



F(d, V ) = ~\C'X - \c'b(Vo + V) 2 , (23) 



where 



_ dCg(d) _ dc b (d + A) 

; °6 — > 



dd 



and Vo can be varied, with V c a fixed property of the 
plates. The force is minimized when 



dF(d,V ) 



dV 



= 



c b v c 



v =v m 



Vm{d) = (25) 



V m {d) 



-Vc 



d 2 + (d + Ay 



eoA 



v; 2 



2 d 2 + (d + Ay 



(27) 
(28) 



Alternatively, in terms of V m (d) (up to V\, sec below), 
the force is 



P ei un toA V 2 (d)[d 2 + (d + A) 2 } 

tresW ~ — ^4 • ( iy J 



Experimentally, V m (d) must include a distance- 
independent offset V\ which arbitrarily depends on the 
sum of contact potentials in the complete circuit be- 
tween the plates. Therefore, the force due to large 
patches considered in this section should be written as 
proportional to (V m (d) + Vi) 2 , instead of simply V 2 n (d). 
where V\ is determined by a fit to experimental data. 
Note that when the two upper capacitor plates are at 
the same distance from the lower capacitor plate, i.e. 
when A = 0, Eq. (|2"T| predicts a V m independent of dis- 
tance, and Eq. (|2T))) predicts a residual electrostatic force 
F°c S (d) oc (V m + Vi) 2 /d 2 . Since no residual force is ex- 
pected in this case, the minimizing potential must be 

Vm = ~Vl 

In order to analyze the sphere-plane geometry, one can 
divide the sphere into infinitesimal planar areas (as done 
in the proximity force approximation), each with a ran- 
dom potential. In this picture, one can think of the two 
upper capacitor plates in Fig. 1 as one of those infinitesi- 
mal parts of the whole spherical surface, and the distance 
A being a local distance A(r) = d+ R(l — cos£?) reflect- 
ing the effect of the curvature of the spherical surface. 
In this case, C' a (d) = — 2ireoR/d, and the denominator of 
Eq. (f2U)) becomes d 2 . Integrating the force on the lower 
planar plate over the whole spherical surface to get the 
net force leads to a further reduction of the power of d in 
the denominator, leaving the sphere-plane residual elec- 
trostatic force proportional to (V m (d) + V\) 2 jd. Again, 
V\ is a fit parameter that represents a sort of surface 
average potential, plus circuit offsets. 
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IV. TOTAL ELECTROSTATIC FORCE 
RESIDUALS 

We are now in a position to compute the total elec- 
trostatic residual force at the minimizing potential. On 
the one hand, the presence of surface patches small with 
respect to the effective area of interaction in the sphere- 
plane geometry leads to an attractive electrostatic force 
given by Eq. (|13[) . which was derived using the proximity 
force approximation (d <C R) to treat the curvature of 
the spherical plate. This component of the force due to 
stochastic potential patches is clearly independent of the 
applied voltage Vq between the plates, and therefore will 
be present even when the applied voltage is set at the 
minimizing potential, Vo = V m (d). As we have seen in 
Section II, the exact dependence of this force on distance 
varies for different models for the statistical properties of 
the two-point correlation functions of the surface poten- 
tials. However, in the limit of patches much smaller that 
the sphere-plane separation (kd <C 1), all models predict 
a dependency of the form V^ ms /d, where V rms denotes the 
rms voltage fluctuations. 

On the other hand, the presence of surface patches 
large with respect to the effective area of interaction, cor- 
responding to a coherent variation of the potential along 
the spherical surface V(r,ip), leads both to a distance- 
dependent minimizing potential V m {d) and to a residual 
component of the electrostatic force F^ s (d) — F(d,Vo = 
V m (d)), even when the external potential Vo is set at 
the minimizing potential. As seen in Section III, this 
force is due to the interplay between surface patches and 
the curvature of the non-planar surface, has the form 
{V m {d) + Vi) 2 /d, and is in addition to the Vb-independent 
force due to small potential patches. 

By adding the two contributions of small and large 
surface patches (Section II and III, respectively), we can 
determine the form of the residual electrostatic force at 
the minimizing potential. By numerical modeling of dif- 
ferent patch sizes and distributions using Eq. (j2"2")l . we 
found that the experimental observations should be de- 
scribed by a relationship of the form 



f£(d) = 7rJk 



(F m (d) + Fi) s 



V 2 



(30) 



where we have specifically assumed the small kd limit for 
the small patches. Both V\ and V ms are constants de- 
termined by fitting to the observed force at large plate 
separations, where the Casimir force is relatively small. 
In flfl ] we have measured the dependency of the min- 
imizing potential with distance, and we have applied 
Eg. ((3TH) to fit the observed residual electrostatic force in 
a Ge sphere-plane Casimir experiment for distances large 
enough (d > 5/xm) to neglect possible contributions due 
to Casimir forces. With only two adjustable parameters 
V\ and Kms, very good fit (% 2 of order unity) between 
the experimental data and Eq. (f3"0")) was possible. We em- 
phasize that Eq. (|30[) was obtained in a heuristic way by 
numerically analyzing the distance-dependent force and 



variation in minimizing potential due to random surface 
potential distributions by use of Eq. (j2"2"]l . The physical 
basis of Eq. ([50)1 is nonetheless quite sound. 



V. DISTANCE CORRECTION DUE TO 
RESIDUAL FORCES AND ITS IMPLICATION 
ON THE INTERPRETATION OF PREVIOUS 
DATA 

In previous work Q, the minimizing potential was as- 
sumed to be independent of distance and the absolute 
separation between the sphere-plane plates in a Casimir 
experiment was determined by fitting a background elec- 
trostatic force (at distances where the Casimir force is 
relatively small and negligible) to a function 



D 



f 



d + d 



of 



(31) 



where Bj and do/ are fitting constants. Unfortunately, 
in [§[ , a possible distance variation of the minimizing po- 
tential was not measured. However, a relatively large 
fixed potential was constantly applied to give a 1/d force 
in the data, and this was used to determine the absolute 
separation. 

One can estimate a systematic error for Q (and pos- 
sibly for other experiments where V m was assumed to be 
constant) by considering a possible distance-dependence 
of the minimizing potential, and taking our Ge measure- 
ment results P3 as "typical" . In the previous section we 
noted that the residual long range force in the Ge exper- 
iment is well described by Eq. (j3"0")) in terms of the mea- 
sured distance-dependent contact potential V m (d). If, 
instead, we had mistakenly assumed V m to be indepen- 
dent of distance and neglected any V rms contribution, the 
same experimental data for the residual long range force 
in fioj is also well described by a power-law fit of the form 
l/d e , with exponent e = 0.72 ±0.2 (as mentioned in foot- 
note [24] in [III). Note that this is the residual force at 
the minimizing potential, so if the applied voltage is not 
set precisely to V m (d), there can be an additional force, 
but we neglect that here. In the following we will as- 
sume that a power-law form for the long-range residual 
force also applies to the Au experiment Q, and we will 
take a nominal value of 0.8 for the exponent to assess 
the possible systematic error introduced in . This esti- 
mate is entirely heuristic, and given that the minimizing 
potential was not measured as a function of distance in 
[8j , it is the best that can be done in a post-analysis and 
represents a reasonable range. 

An error in the distance determination is introduced 
if an improper function 1 j d is used instead of the "true" 
function, which we take to be 1/d 8 . Generally, we might 
expect contributions from both, as 



Bo 



d + d (d + do) ai 



(32) 
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where the first term represents a patch potential force 
and an inexact cancelation of the average contact poten- 
tial (Vi), and the second term represents the force due to 
the variation in V m (d). We assume that do is fixed and 
known; our goal is to determine the distance error in de- 
termining do by fitting to a function of the wrong form. 
We can define an effective x 2 and minimize its value to 
find the offset due to the possibility that the wrong fitting 
function, l/(d+ do) alone, was used in Q. A form that 
allows simple numerical calculation is as follows, substi- 
tuting x = d/do 



dx 



a (1 — a) 

+ (1 + e) ~ (x + l)°-» ~ x + 1 



(33) 

where do/ = do + edo, a parameterizes the relative 
amounts, at d = (or x = 0) of the 1/d and 1/d 08 
forces in the "true" (and assumed known) function, and 
a^max — ^max/do is the largest (dimensionlcss) separa- 
tion to which measurements are taken. When a = the 
"wrong" fitting function is equal to the "true" function, 
and e = 0. For simplicity, we choose to define a as the 
relative contribution of the two forces at do because their 
relative size is distance-dependent. 

In [j| the relative contributions from an applied fixed 
voltage 1/d force, and the possible residual unaccounted 
1/d ' 8 force implied by our recent results [Toj . are in the 
ratio approximately 10/1, so a 0.1, at distances of or- 
der 1 /im. The parameters Bf and e are then determined 
by minimizing x 2 . Although the integral in Eq. (|33[) ap- 
pears as elementary, its evaluation is quite cumbersome. 
Results of numerically minimizing x 2 as a function of a 
and Bf, with x max = 10, show that 



e = 0.65a. 



(34) 



Therefore, for a = 0.1 wp obtain e = 0.065. Applying 
this result to Fig. 4 in [8| shows that the x-axis needs 
to be shifted to the left (toward the origin) by 0.065 /im. 
That is, the distance scale is offset by 0.065 fim. With 
this displacement, the Casimir force F c at the true loca- 
tion is better described by the theoretical result due to 
Bostrom and Sernelius [18j . which can be seen from the 
fractional change in force, 8F C /F C = — 35d/d , which gives 
a 20% effect and brings the measured Casimir force into 
agreement with the predicted result in [l8[ for distances 
around 1 fim; however, this result should be considered 
as preliminary and as a rough estimate of the correction 
magnitude that is possible. Because patch potential ef- 
fects are sample dependent, it is not possible to say con- 
clusively that the effects described here contributed to 
the result, although, not having tested for such possible 
effects, an additional systematic error could be ascribed 
to the result in 

Other possible background forces that deviate from a 
1 jd character will lead to corrections to the true distance 
when a fit to l/(d — do) is performed. For the discussion 
here, we chose a form that is motivated by our experi- 



mental results, and by a theoretical analysis |{|. Alter- 
natively, if a residual force that appears to have a 1/d 
character is removed from Casimir data, a possible back- 
ground force as considered here will persist as a direct 
systematic. 



VI. CONCLUSIONS 

We have derived in a straightforward and heuristic 
manner several important results pertaining to the ex- 
cess electric force between plates that results from ran- 
dom surface patch potentials. These results have been 
cast in terms of the surface autocorrelation function, or 
alternatively the two dimensional spatial Fourier power 
spectrum. Our recent measurement fioj of short-range 
forces with Ge plates in the sphere-plane configuration 
has shown the importance of assessing surface patch po- 
tentials when measuring Casimir force residuals. A re- 
cent Casimir experiment in the plane-plane configuration 
fl3j has also found a large residual (non-Casimir) force 
that is also probably related to electrostatic patch resid- 
ual forces. 

Furthermore, we have shown that long-range surface 
correlations in the sphere-plane geometry can lead to a 
distance dependence of the electrostatic force minimiz- 
ing potential. This effect is due to the dependence of the 
net surface averaging area on the separation between the 
plates and has been described by use of a simple capacitor 
model. The model clearly reproduces the general effects, 
provides an explanation of the origin of the varying min- 
imizing potential, and demonstrates that even when the 
force is minimized, a residual electrostatic force remains. 



The results here should be compared to earlier work [17 1 
where it was assumed that the variation in minimizing 
potential was due to a voltage in series with the plates, 
with that voltage varying with distance by some unknown 
mechanism; this analysis could not describe the observed 
(non-Casimir) force in the germanium measurements de- 
scribed in |l(|. We note that the effect is expected to 
exist whenever there are surface potential patches and 
non-perfectly-parallel surfaces; the model shown in Fig. 
1 requires only the existence of patches with different 
absolute distances between the surfaces. Such distance 
differentials can be created by surface roughness and/or 
lack of parallelism. 

Finally, we have shown that determining the distance 
between a sphere and a plane in a Casimir force exper- 
iment can be subject to systematic effects arising from 
residual electrostatic forces. The magnitude of the error 
in the distance determination is large enough to bring the 
results presented in 0] into agreement with the calcula- 
tion that takes into account properly the low frequency 
permittivity of metals |l8[ . The relevance of the analysis 
presented in this paper to precision Casimir force mea- 
surements and their possible systematic contamination 
should not go unnoticed. We are currently performing 
new measurements using Au coated plates in the appa- 
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ratus used for Ge measurements, and will revisit these 
systematic effects. 
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